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ABSTRACT 

Spectra from the Short Wavelength Spectrometer (SWS) on ISO exhibit ar¬ 
tifacts at 4.5 and 8 /iin. These artifacts appear in spectra from a recent data 
release, OLP 10.0, as spurious broad emission features in the spectra of stars 
earlier than ~F0, such as a CMa. Comparison of absolutely calibrated spec¬ 
tra of standard stars to corresponding spectra from the SWS reveals that these 
artifacts result from an underestimation of the strength of the CO and SiO molec¬ 
ular bands in the spectra of sources used as calibrators by the SWS. Although 
OLP 10.0 was intended to be the final data release, these findings have led to an 
additional release addressing this issue, OLP 10.1, which corrects the artifacts. 

Subject headings: astronomical data bases: miscellaneous — methods: data anal¬ 
ysis — techniques: spectroscopic 


1. Introduction 

The Short Wavelength Spectrometer (SWS) on the Infrared Space Observatory (ISO) 
obtained approximately 1250 spectra covering the full 2.4-45 /irn wavelength range at mod¬ 
erate resolution. We are engaged in an ongoing project to classify these spectra (Kraemer 
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et al. 2001), to reprocess them, and to present them in a publically available database 
on-line (Sloan et al. 2002, in preparation). SWS data are currently available in a partially 
processed form called Auto-Analysis Results (AARs). The AARs have a fairly complex for¬ 
mat (Leech et al. 2001) that requires further processing before the data can be used for 
scientific analysis. Our processing method reduces these data to a single continuous 2.4 to 
45 gm spectrum. 

In assessing the quality of our new reprocessing algorithm, we compared the results 
for several infrared standards observed with the SWS to absolutely calibrated spectra from 
Cohen et al. (1992a, 1992b, 1995, 1996a, 1996b, 2001 (in preparation)). These standards 
are based on synthetic spectra of the A0 dwarf a Lyr and the Al dwarf a CMa, which 
serve as the reference standards in the system. Cohen et al. (1992a) describe the details 
of their method that used high-quality ground-based and airborne photometry to normalize 
the synthetic spectra to measured astronomical fluxes. The synthetic spectra are based on 
the models of Kurucz (1979), with updated opacities and metallicities (Castclli & Kurucz 
1994). 

Secondary standards are added to the system by observing their spectra in conjunction 
with that of the primaries so that atmospheric, telescopic, and instrumental transients can 
be removed. The spectra for the secondary standards are then obtained by dividing the 
observed secondary by the observed reference standard and multiplying by the assumed 
spectrum (the model) for the reference, degraded to match the spectral resolution of the 
instrument used: 
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where the subscript a refers to the reference standard, and the subscript b refers to the 
secondary. This is the standard method used by spectroscopists at ground-based telescopes 
to calibrate a program source by ratioing its spectrum to that of a standard star, preferably 
secured at an airmass matching that of the target. In effect, Equation 1 transfers the quality 
of the synthetic A star model from the reference to the new standard, whatever its spectral 
type. 

Cohen et ah (1992b, 1995, 1996a, 1996b, 2001) have applied this method to create 
absolutely calibrated composite spectra for 13 infrared standards, using spectra from ground- 
based telescopes, the Kuiper Airborne Observatory (KAO) and the Infrared Astronomical 
Satellite (IRAS) (further details of the process can be found in Cohen et al. 1992a,b). 
Cohen et ah (1996a) added a 1 Cen (G2 V) as an alternative reference standard for the 
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southern hemisphere to give a total of 3 reference standards as well as the 13 secondaries. 
Most of the secondary standards are giants with spectral classes later than KO and were 
chosen for their intrinsic brightness; sources later than approximately M3 are avoided due 
to the possibility of their variability (Eyer & Grenon 1997). 

The preferred infrared reference standard is Sirius (a CMa), due to its brightness and 
its dust-free spectrum beyond 20 /im. Figure 1 compares our derived SWS spectrum of this 
source to the Kurucz model presented by Cohen et al. (1992a). Deviations between the two 
occur in the vicinity of 4.5 //rri and again at 8.0 pm. We propose an explanation for the 
origin of this discrepancy in this Letter and discuss the implications for the calibration of 
the SWS database and the impact on calibration of future infrared missions. 


2. Data Format and Analysis 

Data from the SWS are publically available in AAR format. This format corrects the 
spectra, as far as possible, for the relative spectral response of the detectors, differences in 
gains between individual detectors and a variety of other issues, as described by Leech et al. 
(2001). The hies are organized in 12 spectral segments, ranging from 1 to 4 for each of the 
4 different detector bands. Each spectral segment contains interleaved spectra from the 12 
individual detectors, scanned in the direction of both increasing and decreasing wavelength, 
giving a total of 24 separate spectra for each spectral segment. Thus the processing of the 
full wavelength range from an AAR hie requires the combination of 288 individual spectra 
into one. 

We start with data from the Off-line Processing (OLP) pipeline 10.0, available in June 
2001, which was intended to be the final release. This version includes the latest attempts 
to remove the memory effect from the Si:As data (4-12 pm) (Kester et al. 2001). Although 
these attempts have significantly improved the data, residual effects remain, and one can see 
that the spectra presented here tend to diverge from the expected result as the data approach 
12 pm. However, we have high confidence in the quality of the data at shorter wavelengths, 
where the flux levels from the sources are much stronger and where the discrepancies between 
the SWS data and the reference spectra occur. 

Sloan et al. (2002) will present details of the algorithm used to combine the 288 in¬ 
dividual spectra into one. We have tested for the possibility that the algorithm might be 
responsible for the artifacts seen in Figure 1 by processing the AAR data manually with 
the Infrared Spectral Analysis Package (ISAP) following the standard procedures. The two 
methods produce similar results, which indicates that the problem is intrinsic to the spec- 
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tram released by ISO and not to a specific method of processing. The artifact peaking at 
4.5 /ini also appears in Figures 6 and 7 of Decin (2001), although the wavelength stretch 
makes it more difficult to notice. 


3. Discussion 

The divergence between the SWS spectrum of a CMa and the Kurucz model is char¬ 
acterized by two broad “emission” features, at ~4-5 //rn and ~8 /in i. (The rise at longer 
wavelengths is caused by the memory effect mentioned in Sec. 2.) The feature at 8 //rn is 
reminiscent of the 8 //in emission artifact found in the data from the Low Resolution Spec¬ 
trometer (LRS) aboard IRAS by Cohen et al. (1992b). The LRS database was originally 
calibrated assuming that a Tau could be characterized as a 10 4 K blackbody. However, since 
the spectrum of a Tau (a K5 III) contains significant absorption from the SiO fundamental 
band at 8 //rn, this error propagated to the entire LRS database, producing an apparent 
emission feature in the spectra of a CMa, a Lyr, and other early-type stars. 

Fifteen stars have served as spectral calibrators for the SWS (Schaeidt et al. 1996; Leech 
et al. 2001; Shipman et al. 2001). Ten of them are K and M giants, which show strong 
molecular absorption bands in their spectra. The dominant bands peak at wavelengths of 
~2.4 /tm (CO overtone), ~4.1 /tm (SiO overtone), ~4.5 //rn (CO fundamental, blended with 
the weaker SiO overtone), and ~8.0 /tm (SiO fundamental). For the assumed spectrum of 
each of these stars in OLP 10.0, the SWS instrument team chose to use synthetic models 
of the stars from 2 to 12 /tm. For longer wavelengths, a composite spectrum from Cohen et 
al. where available or a template spectrum from Cohen et al. (1999) where no composite 
existed were used. The two sets of data were spliced together at 12 /tm (Schaeidt et al. 
1996; Leech et al. 2001; Shipman et al. 2001). Previous pipeline releases used the Cohen 
composites and templates for the entire wavelength range. 

The K5 giant 7 Dra is a typical example of an SWS calibrator. Figure 2 compares 
the synthetic spectrum of this source (Decin 2000, 2001) to the observationally calibrated 
composite spectrum (Cohen et al. 1996b). The differences between the two spectra (the 
lower panel) are similar to the artifacts seen in Figure 1. Decin (2001) compares the 
composite spectra of several other SWS calibrators to the synthetic spectra (in her Figure 
9), and most of them show similar deviations 6 . It is highly unlikely that an instrumental 


(> Decin (2001) notes the discrepancy between the model and the composite spectrum of 7 Dra, but 
incorrectly attributes it to the substitution of data of a Tau for 7 Dra in this wavelength regime. In the 
header to the file containing the composite spectral data for 7 Dra, Cohen et al. (1996b) state that they 
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Fig. 1.— A comparison of SWS data for a CMa (dotted line) with the Kurucz model used 
by Cohen et ah (1992a) (solid line). The spectra have been normalized between 6 and 7 /im 
and are plotted in units of A 4 F\, which clearly shows deviations from the reference spectrum 
at the CO and SiO bands. The bottom panel plots the difference between the two as a 
percentage of the flux density of the Kurucz model. 
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effect (such as the memory effect) would manifest itself precisely at the wavelengths of these 
photospheric absorption bands. 

The coincidence of the emission artifacts apparent in Figure 1 with the molecular bands 
in K and M stars and the similarity of the artifacts to the deviations between the model 
and composite spectrum in Figure 2 lead us to hypothesize that the artifacts result from an 
underestimate of the depth of these bands in the synthetic spectra. 

Figure 3 compares the data from the SWS to both the model and composite spectrum 
of 7 Dra. Although not perfect, the SWS spectrum, as calibrated in OLP 10 . 0 , matches the 
model of 7 Dra significantly better than the Cohen composite. Additional spectra of 7 Dra 
from the SWS database, taken at different times during the mission and at different spectral 
resolutions, as well as spectra of other cool calibrators, produce similar results. 

If the synthetic spectra used to calibrate the SWS generally underestimate the strength 
of the molecular bands, then this miscalibration should propagate to the entire SWS database. 
As Figure 4 shows, the artifacts at 4.5 and 8 /mi are readily seen in bright spectra from early- 
type stars, including the reference standard a Lyr, which should show neither emission nor 
absorption from molecular bands. In the spectrum of the average giant with no compos¬ 
ite spectrum for comparison, the effect of the miscalibration will be more subtle and more 
difficult to recognize, since it will only reduce the apparent depth of the absorption bands. 

In Figure 3 the difference between the SWS data and the composite spectrum for 7 Dra 
exhibits the largest excess between ~4.0 and 4.7 /mi. While an excess clearly exists at these 
wavelengths in the spectra of the hot stars in Figure 4, the confidence level of the magnitude 
of this excess shown in Figure 3 is not high. Because the atmosphere is almost opaque in the 
core of the CO 2 band at 4.3 /jrn, even at aircraft altitudes, the uncertainty in the spectral 
composite is typically ~ 6 % from 4.22 to 4.58 /jrn and peaks at 10% at 4.22 /mi (Cohen et al. 
1996b). Roughly half of the composites, including 7 Dra, have a data gap at ~4.2-4.4 /mi 
due to the atmospheric band. The spectral artifact is, however, much broader than the data 
gap and larger than the uncertainties can account for. The uncertainties in the composite 
spectrum at the other feature are < 2 %. 

To reduce the noise from a single comparison, we calculated the weighted average of the 
ratios of all of the spectra common between the sources used to spectroscopically calibrate 
the SWS and the composite and reference spectra from Cohen et al. This weighted mean, 
displayed in Figure 5, provides an estimate for the modifications needed in the OLP 10.0 
data. It also provides an assessment of the reliability of the synthetic spectra relative to the 


calibrated this spectral region using spectral ratios of 7 Dra to a Boo obtained on the KAO. 
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absolutely calibrated composite spectra. 

In light of our findings, the SWS team has revised their calibration strategy and the ISO 
Data Center has agreed to an additional data release, OLP 10.1. By basing the new spectral 
calibration only on stars earlier than K, the worst of the discrepancy between the synthetic 
spectra and the Cohen spectra should be mitigated. The Decin models agree with the Cohen 
composites, templates, or Kurucz models (as appropriate) for the five stars (two As, one F, 
and two Gs) which meet this criterion (R. Shipman 2001, private communication). 


4. Conclusions 

We have identified spectral artifacts in the vicinity of the CO and SiO bands of late- 
type standard stars that appear in all OLP 10.0 spectra. These artifacts constitute a known 
systematic bias over a fairly large spectral range, and removing them is relatively straight¬ 
forward. While OLP 10.0 was to have been the final version, an additional release, OLP 
10 .1, has been made to correct the database based on our findings. 

The SWS calibrations between 2 and 12 /irn were predominantly based on stellar atmo¬ 
spheric models of cool stars. These models have been used with great success in analysis of 
physical properties of cool stars (e.g. Decin et ah 2000). Model spectra of A stars, which 
contain only atomic lines and no molecular bands, appear to be well founded. Synthetic 
spectra can also achieve higher spectral resolution than that afforded by the composites. 
However, the artifacts in the SWS database that mirror the dominant molecular absorption 
features in the 4-10 /irn range indicate that synthetic spectra of cool giants require further 
progress before they can be used for definitive spectral calibration. 

The Infrared Spectrograph (IRS) on the upcoming Space Infrared Telescope Facility 
(SIRTF) faces issues similar to those encountered by the SWS team, with the added com¬ 
plication that all of the commonly used spectral standards will be far too bright for use 
on the IRS. Until the commission of the Stratospheric Observatory for Infrared Astronomy 
(SOFIA), it will prove difficult to observationally calibrate standards faint enough for use by 
the IRS. Until this happens, the IRS will need to rely on models for spectral calibration. If 
models of later stars are required, improvements in the production of synthetic spectra will 
also be necessary. 


We wish to thank Russ Shipman for many useful discussions regarding the calibration 
process. Thijs de Graauw provided support for this work through the USAF/SRON MOLI. 
One author (KEK) was supported by the National Research Council via a Research Associ- 
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Fig. 2.— A comparison between the model of 7 Dra used to calibrate the SWS data (thin 
line ) and the composite spectrum calibrated by Cohen et ah (1996b) (thick line). The 
spectra have been normalized between 6 and 7 pm and the model resampled to match the 
composite. The data for the model have been adapted from Fig. 8 and 9 in Decin (2001). 
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Fig. 3.— Comparing SWS data for 7 Dra (dotted line ) with calibrator spectra. The upper 
panel is as in Fig. 2, but with a spectrum from the SWS superimposed (dotted line). The 
lower panel compares the difference between the SWS data and the model (thin line) to the 
difference between the SWS data and the composite from Cohen et al. (1996b) (thick line). 
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Fig. 4.— Several SWS spectra of bright early-type naked stars. Most have been multiplied 
by a constant (given in parentheses). The bottom spectrum is the Kurucz model of a Lyr 
used by Cohen et al. (1992a). 
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Fig. 5.— The weighted mean of the deviations between the SWS spectra of the standard 
stars in common with the composites and models of Cohen et al. Both the SWS data and 
the composites were resampled to 0.1 /jrn spacing and placed on the same wavelength grid. 
Error bars are one sigma. 
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